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Abstract 

We construct a supersymmetric 50(10) xU(1)f model of the Yukawa interactions at 
the grand unification scale from knowledge of a phenomenological set of mass matrices 
obtained by a previous bottom-up approach. The U(1)f family symmetry determines 
the textures for the Majorana and generic Dirac mass matrices, while the 50(10) 
symmetry relates each particular element of the up, down, neutrino and charged lepton 
Dirac matrices. The dominant second and third family contributions in the Dirac 
sector are renormalizable, while the remaining contributions to the Dirac mass matrices 
are of higher order, restricted by the U(l)p family symmetry to a small set of tree 
diagrams, and mainly complex-symmetric. The tree diagrams for the Majorana mass 
matrix are all non-renormalizable and of progressively higher-order, leading to a nearly 
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geometrical structure. Pairs of 1 , 45,10 and 126 Higgs representations enter with those 
having large vacuum expectation values breaking the symmetry down to SU( 3) c x 
SU( 2)l x U(1)y near the grand unification scale. In terms of 12 parameters expressed 
as the Yukawa couplings times vacuum expectation values for the Higgs representations 
employed, a realistic set of 15 quark and lepton masses (including those for the 3 heavy 
righthanded Majorana neutrinos)and 8 mixing parameters emerges for the neutrino 
scenario involving the non-adiabatic conversion of solar neutrinos and the depletion of 
atmospheric muon-neutrinos through oscillations into tau-neutrinos. 

PACS numbers: 12.15.Ff, 12.60Jv 
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I. INTRODUCTION 

The standard model (SM) of strong and electroweak interactions, while providing excel¬ 
lent agreement with experiment todate, is known to be woefully inadequate to explain the 
mass spectrum and mixings of the three families of quarks and leptons. One needs to go be¬ 
yond the standard model in order to relate the independent Yukawa couplings to each other. 
Of the various possibilities, supersymmetric grand unified theories and superstring theories 
seem to hold the most promise for successfully incorporating the Yukawa interactions in a 
more satisfactory fashion. In this paper we shall restrict our attention to supersymmetric 
SO (10) grand unification, which has been shown [1] to unify the gauge couplings successfully 
at a scale of A sgut rs _/ 10 16 GeV. 

It is a generally held opinion that knowledge of the mass matrices in the weak flavor 
basis can provide insights into the dynamical mass-generating mechanism. [2] This follows 
from the fact that the mass eigenvalues are obtained by diagonalization of the mass matrices, 
while the mixing matrices in the mass eigenbasis can be constructed from knowledge of the 
diagonalizing matrices connecting the two bases. By starting from the correct mass matrices, 
one should then be able to deduce the observed quark and lepton masses and mixings after 
evolving the results down to the present “low energy” scales. 

Generally two procedures are at one’s disposal for the identification of the “correct” 
mass matrices. One can attempt to postulate a particular structure or “texture” for the 
mass matrices based on some well-defined and presumably simple theoretical concepts such 
as the unification group and/or the number of texture zeros present.[3] This procedure has 
been employed by most researchers in the past twenty years. Alternatively, one can make 
use of the known low energy mass and mixing data, supplemented by reasonable guesses 
for data which is not yet well determined, in order to extract mass matrices within some 
framework at the unification scale which yield the low energy data in question. Of special 
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interest are neutrino scenarios incorporating the Mikheyev - Smirnov - Wolfenstcin (MSW) [4] 
nonadiabatic resonant conversion interpretation of the depletion of solar electron-neutrinos 
[5] and either the observed depletion of atmospheric muon-neutrinos through oscillations [6] 
or neutrinos of satisfactory mass to contribute to the hot component of mixed dark matter 
[7], for example. 

In a series of papers [8] the authors have demonstrated the latter “bottom-up” approach 
by making use of Sylvester’s theorem [9] to construct mass matrices from the low energy 
masses and mixing matrices evolved to the unification scale. In doing so, we have attempted 
to look for simplicity of the mass matrices in the 50(10) framework while varying the quark 
and lepton weak bases. Such simplicity was found for the MSW solar and atmospheric 
neutrino depletions in the bases where the up quark and Dirac neutrino mass matrices are 
real and diagonal, while the down quark and charged lepton matrices are in general complex 
symmetric. The right-handed Majorana neutrino matrix exhibits a simple nearly geometrical 
texture. 

From the phenomenological mass matrices constructed, we have attempted to derive mass 
matrices of similar textures from some well-defined family symmetry. In particular, we find 
within an 50(10) x U{1)f symmetry framework that we can reproduce all the known and 
assumed-known low energy mass and mixing data for the quarks and leptons in terms of 
products of Yukawa couplings and Higgs vacuum expectation values (VEVs). The U(l)p 
symmetry controls the textures for the generic Dirac and Majorana mass matrices, while 
50(10) relates particular elements of the up, down, Dirac neutrino and charged lepton mass 
matrices to each other. 

In this paper we shall present all the details for this model construction which were 
summarized in a short letter submitted elsewhere [10]. Section II summarizes the bottom- 
up procedure and the phenomenological mass matrices obtained for the neutrino scenario 
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preferred. The U(1)f family symmetry is introduced and applied in the Dimopoulos tree- 
diagram approach [11] in Sect. Ill for the contributions to the mass matrices. In Sect. IV the 
diagramatic contributions to the mass matrix elements are explicitly given with quantitative 
results presented in Sect V. Our work is summarized in Sect. VI. 

II. PHENOMENOLOGICAL MATRICES from a BOTTOM-UP APPROACH 

We begin by presenting the low scale input and procedure by which we were able to 
construct a relatively simple S'O(IO) set of phenomenological mass matrices at the SUSY 
GUT scale as spelled out in detail in Ref. [8] The relevant framework is assumed to be that 
of SUSY £0(10) grand unification at a scale of A sgut = 1-2 x 10 16 GeV with supersymmetry 
breaking occuring at a scale of 180 GeV, in order that we can use the analytical one-loop 
evolution formulas and results given by Naculich [12]. 

For the low scale quark data, we assumed the following set of quark masses and Cabbibo 
- Kobayashi - Maskawa (CKM) mixing matrix [13] 


m„(lGeV) 

= 5.1 MeV, 

m d (lGeV) = 8.9 MeV 


m c {m c ) 

= 1.27 GeV, 

m s (lGeV) = 175 MeV 

(2.1a) 

m t {m t ) 

= 150 GeV, 

m b (m b ) ~ 4.25 GeV 



/ 0.9753 

0.2210 0.0031e _,155 ° \ 



VcKM = 

-0.2206 

0.9744 0.043 


(2.16) 


0.011 

-0.041 0.999 




^ —0.001* 

) 



The light quark masses were chosen to be the central values given by Gasser and Leutwyler 

[14], while the heavy physical top mass was 

set equal to 160 GeV 

prior 

to its discovery 

yielding a running mass of 150 GeV. We assumed a value of 0.043 for V cb , which is now 


thought to be closer to 0.040, and applied strict unitarity to determine V ub , V t d and V ts . 
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The greatest AO (10) simplicity was obtained for the neutrino scenario incorporating the 
observed depletion of solar neutrinos [5] through the nonadiabatic MSW [4] matter conversion 
of electron-neutrinos into muon-neutrinos in the interior of the sun and the depletion of 
atmospheric muon-neutrinos through oscillation into tau-neutrinos observed now by several 
deep mine collaborations [6]. The central values deduced for these mixing plane results are 


Sml 2 

~5x 10“ 6 eV 2 , 

sin 2 2d V2 

~ 0.008 

5m 2 23 

~ 1 x IQ’ 2 eV 2 , 

sin 2 2023 

~ 0.9 


We took for the lepton input 


( 2 . 2 ) 


m Ue = 0.5 x 10 6 eV, m e = 0.511 MeV 

m Vfl = 0.224 x 10“ 2 eV, = 105.3 MeV (2.3a) 

m VT = 0.105 eV, m T = 1.777 GeV 

and 

/ 0.9990 0.0447 0.0076e- i155 

T/ -0.0363 0.8170 0.575 

* lept 

0.026 -0.570 0.818 

V —0.007i 

These masses and mixing matrix data were evolved to the SUSY GUT scale by using 
formulas given by Nacnlich [12] as spelled out in detail in Ref. 8. We could then reconstruct 
complex-symmetric mass matrices at the SUSY GUT scale by making use of Sylvester’s 
theorem [9] as illustrated by Kusenko [15] for the quark sector. The construction is not 
unique, for one is free to change the quark and lepton weak bases by letting two parameters, 
x q and xg, vary independently over their support regions, 0 < x < 1. For x q (xg) = 0, the 
up quark (Dirac neutrino) mass matrix is diagonal; while for x q (xg) = 1, the down quark 
(charged lepton) mass matrix is diagonal. One is also free to vary the signs of the mass 
eigenvalues. 



By varying the signs of the mass eigenvalues and the two parameters x q and xg, we 
then searched for a simple S'O(IO) structure for the mass matrices. The greatest simplicity 
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occurred with x q = 0 and xt, = 0.93 corresponding to diagonal up quark and Dirac neutrino 
mass matrices and leading to 


M u ~ M N ° irac ~ diag{ 126; 126; 10) 



/ 10', 126 

10', 126' 

10'\ 

m d ~ m e ~ 

10', 126' 

126 

10' 


V 10' 

10' 

10/ 


(2.4a) 


(2.46) 


with M E , and M|[ anomalously small and only the 13 and 31 elements complex. En¬ 
tries in the matrices stand for the Higgs representations contributing to those elements, 
which we elaborate upon in the next Section. We have assumed complete unification for the 
Yukawa couplings of the third families of quarks and leptons and that vacuum expectation 
values (VEVs) develop only for the symmetric representations 10 and 126. The 10’s con¬ 
tribute equally to (M u , M D ) and (. M Noirac , M E ), while the 126’s weight (M u , M D ) and 
(.M x,niac \ M e ) in the ratio of 1 : -3. The Majorana neutrino mass matrix M R , determined 
from the seesaw formula [16] with use of M N ° irac and the reconstructed light neutrino mass 
matrix, exhibits a nearly geometrical structure [17] given by 



/ F 

-Vfe 

Vfc \ 

m r ~ 

-Vfe 

e 

-Vec 


{ VFC 

-Vec 

c ) 


(2.4c) 


where E ~ | \JFC with all elements relatively real [18]. It can not be purely geometrical, 
however, since the singular rank-1 matrix can not be inverted as required by the seesaw 
formula, M Neff ~ —M NDirac (M R )~ 1 M N Dirac. 


III. U(1) F FAMILY SYMMETRY and RESULTING TREE DIAGRAMS 

The challenge is now to introduce a family symmetry which will enable us to derive the 
mass matrix textures derived above phenomenologically from our bottom-up approach. For 
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this purpose, we propose to use a U(1)f family symmetry [19], where we leave open for the 
time being whether the symmetry is global or local in which case it can be gauged. Before 
proceeding with this, we review briefly the elements of the AO(10) symmetry group which 
play important roles in our model construction. 

In the SO(IO) framework, each family of left-handed quarks, leptons, conjugate quarks 
and conjugate leptons belongs to a 16 dimensional representation. It is convenient to rep¬ 
resent a given flavor (and color) member of the ith family and its conjugate by the two 
components v k,x = (tpn, {'fp c )iL )• In the corresponding three-family basis ordered as follows, 
V L — {ipiLi ip2L-> ip3L, Li {' l P°)2L, {' 4 ’°)?,l} the contributions to the np or down quark, neu¬ 
trino or charged lepton, mass terms in the Yukawa Lagrangian are then given by 

C = VlC- 1 MV L + h.c. (3.1a) 


where the 6x6 matrix can be written in terms of 3 x 3 submatrices 

M l Mr 


M = M T = 

with the individual contributions referring to 


M t 

Dirac 



M L : (AlVC-^jl 
MDirac ■ (V’»l) T C' _ 1 (V’ C )jL = 'IpjR^iL 

M Dirac : {^ C )iL T C-^jL = 'IpiR'lpjL 

M R : {^ c )J L C- l {ip c ) jL = VkR {$ c ) jL 


( 3 - 16 ) 


(3.1c) 


Here the diagonal block entries appear only for neutrinos with M L the left-handed Majorana 
neutrino mass matrix which we take to vanish, while M R is the right-handed Majorana 
neutrino mass matrix which receives large contributions near the SUSY GUT scale. 


By construction the 6x6 matrix A4 is complex symmetric, but the Dirac mass submatrix 
is not necessarily complex symmetric. We shall assume that the dominant contributions are 
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complex symmetric and that any departures from this form arise from small higher-order 
corrections. Recall that the 50(10) product rules read 


16 x 16 

= 10, + 120 a + 126, 

(3.2a) 

16 x 16 

= 1 + 45 + 210 

(3.26) 


Hence we shall assume that only the symmetric Higgs representations 10 and 126 develop low 
scale VEVs, while the antisymmetric 120 does not. In terms of the SU( 5) decompositions, 
we have 

10 —> 5 + 5, 126 —>50 + 45 + 15 + 10 + 5 + 1 (3.3a) 

The up-type quarks and Dirac neutrinos then can receive contributions from the neutral 
members of 10(5) and 126(5), the down-type quarks and charged leptons from those of 
10(5) and 126(45), and the heavy right-handed Majorana neutrinos from those of 126(1). 
We shall later assume the Higgs representations 1 and 45 play a role in the higher-order 
corrections, where the 45(1) and 45(24) develop VEVs according to the decomposition 

45 ^24+ 10+10 + 1 (3.36) 

Returning to the phenomenological mass matrices obtained in Section II, we use the 
textures given in (2.4a,b,c) as our starting point for the construction of an 50(10) x U(l)p 
model of the Yukawa interactions. We find it useful to introduce a generic Dirac matrix, 
M Dirac, to go along with the one Majorana matrix, M R . The U(1)f family symmetry will 
then determine the textures for Moi rac and M R , while the 50(10) symmetry will relate the 
corresponding matrix elements of the four Dirac matrices M u , M D , M N and M E to each 
other. 

Simplicity of the 50(10) structure requires that just one Higgs 10 representation con¬ 
tributes to the (M Dirac ) 3 3 element (hereafter labeled D33). Since a 10 contributes equally to 
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the 33 elements of all four Dirac matrices, this implies that we assume complete unification 
of the Yukawa couplings at the unification scale: m T = m b = rri t / tan/ho, where tan R 10 is 
equal to the ratio of the up quark to the down quark VEVs in the 10, i.e., 

frit = 9io(v/V 2 ) sin (ho = 9ioV u 

m b = rh T = g 10 {v/s/2) cos(3 10 = g w v d (3.4a) 

tan (3 W =v u (5)/v d (5) 

in terms of the SU( 5) decomposition of 50(10) with v = 246 GeV. The same 10 can not 
contribute to D23 = D32, for the diagonal nature of M u and M N requires the presence of 
another 10' with 

tan Pi 0 i = v' u (5')/v' d (5') = 0 (3.46) 

Likewise we assume a pure 126 contribution to D22 with 

tan /%26 = w u (5)/w d (45) (3.4c) 

In other words, for simplicity we have taken the 2-3 sector of Movrac to be renormalizable 
with two 10’s and one 126 developing low scale VEVs. We illustrate the renormalizable 
3-point tree diagrams in Fig. la. 

We now assign U(l)p charges as follows to the three families (in order of appearance) 
and to the three Higgs representations introduced which generate low scale VEVs with the 
numerical values to be determined later: 

16“, 16f, 16j, 10 a , 10' 6 , 126 c (3.5a) 

Conservation of U(1)f charges then requires 2cc + a = 0, a + (3 + 6 = 0 and 2/3 + c = 0 as 
seen from the diagrams in Fig. la. 

We assume the rest of the MDirac elements arise from higher-order tree diagrams as first 
suggested by Dimopoulos [11] twelve years ago. The point is that not only does SUSY 
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control the running of the Yukawa couplings between the SUSY GUT scale and the weak 
scale where it is assumed to be softly broken, but it also allows one to assume that only 
simple tree diagrammatic contributions to the mass matrices need be considered as a result 
of the non-renormalization theorem [20] applied to loop diagrams. While the low-scale VEVs 
introduced act only once in each diagram, other GUT scale VEVs arising from 1 and 45 
Higgs representations can connect superheavy GUT scale 16 fermions and their conjugate 
16 mirrors to each other and to the three light 16 families. The superheavy fermions and 
their mirrors pair off and get masses slightly higher than the SUSY GUT breaking scale, so 
their propagators in the higher-order tree diagrams lead to non-renormalizable contributions 
scaled by their masses. 

For each 45 Higgs representation, as noted earlier in terms of the SU(5) decomposition 
given in (3.3b), VEVs can develop in the orthogonal directions 

< 45 x >~ 45(1), < 45 y >~ 45(24) (3.6a) 

or in any non-orthogonal directions. One such direction of interest corresponds to the hy¬ 
percharge direction for flipped [21] SU( 5) x U( 1) as clarified in Table I: 

< 45 z >= ^ < 45 Y > < 45 r > (3.66) 

5 5 

While the < 45 y > VEV breaks AO (10) -»• SU( 5), the < 45 z > VEV breaks AO (10) -> 
flipped SU( 5). Alternatively, if the AO (10) —■> SU( 5) breaking occurs earlier by some other 
VEV such as < 126 / > as required later for the Majorana sector, the combined action of 
< 45 y > and < 45^ > will result in the breaking of SU( 5) —>• SU( 3) c x SU(2) L x C/(l)y. 

Since the D13 and D23 elements in (2.4a,b) have the same 10' transformation property, 
this suggests that we introduce a 45^ Higgs held and construct an explicitly complex- 
symmetric dimension-6 tree diagram as shown in Fig. lb, for which U(l)p charge conserva¬ 
tion requires a + 7 + 6 + 2e = 0. We shall later give the four Dirac mass matrix contributions 
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derived from D13 by use of Table I which confirms that D13 and D23 do have the same 10' 
transformation property, i.e., the contributions to M u and M N vanish while those to M D 
and M E are non-zero and equal. 

The D12 element, on the other hand, appears to arise from a linear combination of 10' and 
new 126 VEV contributions for which (M E ) 12 <C (M D ) 12 . Rather than introduce another 
new renormalizable diagram, we can make use of the fact that a 45^ Higgs develops a VEV 
which vanishes for the charged lepton D12 diagram as seen from Table I. We then introduce a 
new 45| Higgs field and construct the complex-symmetric dimension -6 tree diagram shown 
in Fig. lb. Note that detailed study showed that to reduce the number of contributing 
diagrams the 10' Higgs line should leave the diagram, or equivalently, the 10'* line should 
enter the diagram, so U(l)p charge conservation requires f3 + 7 — b + 2h = 0. 

The Dll element is dimension -8 or higher, and we leave it unspecified. The complex- 
symmetric leading-order Yukawa diagrams which we wish to generate are then neatly sum¬ 
marized by the ordering of the Higgs fields where all external lines enter the diagrams: 

.D33 : 16 3 -10 - 16 3 
D23 : 16 2 -10' - 16 3 
D32 : 16 3 -10' - 16 2 
D22 : 16 2 —126 — 16 2 

(3.7a) 

D13 : I61 —45 y -10' - 45 x -16 3 
D31 : 16 3 -45 a - -10' - 45 x -I61 
£>12 : I61 —45 z -10'* - 45 z -16 2 
£>21 : 16 2 -45 z -10'* - 45 z -I61 

In order to obtain a different set of diagrams and hence a different texture for the Majo- 
rana matrix, we begin the M33 contribution with a dimension -6 diagram shown in Fig. lc 
by including a new 126 ,J Higgs which develops a VEV at the GUT scale in the SU( 5) singlet 
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direction, along with a pair of l 9 Higgs fields. Here 2a + d + 2g = 0. The nearly geometric 
structure [8] for M R can then be generated by appending more Higgs fields to each diagram. 
For M23 we introduce another 1'^ Higgs field to construct a diagram with one 126 / , one 45 
one IF and two I s fields with charge conservation demanding a + /3 + d + 2g + e + f = 0. 
The new V held is needed in order to scale properly the Majorana matrix elements rela¬ 
tive to each other. The remaining leading-order diagrams of the complex-symmetric Majo¬ 
rana mass matrix follow by appending more 45^, 45| and 1'^ Higgs lines. The pattern is 
made clear from the charge conservation equations: 2/3 + d + 2g + 2e + 2/ = 0 for M22, 
a + 7 + d + 2g + e + h + 2/ = 0 for M13, (3 + 7 + d + 2g + 2e + h + 3/ = 0 for M12, and 
2y + d + 2g + 2e + 2h + 4/ = 0 for Mil. 

In summary, the following Higgs representations have been introduced in addition to 
those in (3.5a): 

l26 ,d , 45^, 45|, l 9 , l' f (3.5 b) 

all of which generate massive VEVs near the GUT scale. In order to obtain CP-violation 
in the quark and lepton mixing matrices, we allow the VEVs for 45x, 45^, 1 and V to be 
complex, but the VEVs associated with the 10, 10' , 126 and 126 representations can be 
taken to be real without loss of generality as seen from our bottom-up results. Clearly, many 
permutations of the Higgs fields are possible in the higher-order diagrams. 

At this point a computer search was carried out to generate U(l)p charge assignments 
leading to the fewest additional diagrams allowed by charge conservation. An especially 
interesting charge assignment stood out for which 


a = 9, (3 — —1, 7 = —8 

a = -18, b = - 8 , c = 2, d = -22, e = 3.5, / = 6.5, g = 2.0, h = 0.5 


(3.8a) 


One should note that since a + /3 + 7 = 0, the [S'O(IO )] 2 x U(1)f triangle anomaly vanishes, 
whereas the [U( 1)f ] 3 anomaly does not. Simplicity then suggests that the U (1 )f family sym- 
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metry group can be global with a familon being generated upon its breaking. Alternatively, 
the U{l)p group can be local and gauged if the [U(1)f] 3 anomaly is canceled by the addition 
of several S'O(10)-singlct fermions with appropriate U( 1) charges, or perhaps better still, by 
the Green-Schwarz mechanism [22] provided the model can be derived from string theory. 
We intend to study this point at greater length elsewhere and do not commit ourselves here 
to either possibility. 

With the above charge assignments we can further greatly limit the number of permuta¬ 
tions and eliminate other unwanted diagrams by restricting the U(l)p charges appearing on 
the superheavy internal fermion lines. With the following minimum set of allowed charges 
for the left-handed superheavy fermions Fp and their mirror partners 

F l : -0.5, 1.0, 2.0, 4.0, 4.5, -4.5, -7.5, 11.0, 12.5 

(3.8 b) 

F c l : 0.5, -1.0, -2.0, -4.0, -4.5, 4.5, 7.5, -11.0, -12.5 

as determined by another computer program, we recover just the leading-order diagrams 
listed in (3.7a) for the generic Dirac mass matrix together with the following uniquely- 
ordered diagrams for the complex-symmetric Majorana mass matrix 

M33 : 16:3 -1 - 126' - 1 - 16 3 

M23 : 16 2 -1 - 45 x -1' - 126' - 1 - 16 3 

M32 : 16 3 -1 - IM' - 1' - 45 x -1 - 16 2 

M22 : 16 2 -1 - 45 x -1' - 126' - 1' - 45. Y -1 - 16 2 

M13 : 16i —45 y -1' - 1 - 45 z -1' - IM' - 1 - 16 3 (3.76) 

M31 : 16 3 -1 - 126' - 1' - 45 z -1 - 1' - 45 x -16, 

M12 : 16! —45 y -1' - 1 - 45 z - 1 ' - IM' - 1 ' - 45 x -1 - 16 2 

M21 : 16 2 -1 - 45jy -1' - 126' - 1' - 45 z -1 - 1' - 45 x -16i 

Mil : 16i -45 x -1' - 1 - 45 z -1' - IM' - 1' - 45 z -1 - 1' - 45 A —16i 

Several other higher-order diagrams are allowed by the U(l)p charges given in (3.8a,b) and 
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appear for Dll, D22, M23 and M32 with the Higgs fields ordered as follows: 

Dll : 16 1 —45 x -1' - 1 - 10' - 1 - V - 45. Y -16i 

D22 : 16 2 -45 z -10'* - 1'* - I 69 , 16 2 -1'* - 10'* - 45 z -16 2 

(3.7c) 

M 23 : 16 2 -45^ -1' - 1 - 45 z -1' - 126 - 1 - 16 3 

M32 : 16 3 -1 - 126' - V - 45 z -1 - V - 45^ -16 2 

These corrections to M23 and M32 ensure that M R is rank 3 and nonsingular, so that the 
seesaw formula [16] can be applied. Up to this point the contributions are all complex- 
symmetric. 

Additional correction terms of higher order which need not be complex-symmetric can be 
generated for the Dirac and Majorana matrix elements, if one allows additional superheavy 
fermion pairs with new U(l)p charges. Such a subset which does not destroy the pattern 
constructed above, but helps to improve the numerical results for the lepton masses and 
mixings, consists of the following: 


F l : 1.5, -6.0, -6.5 

F c l : -1.5, 6.0, 6.5 


(3.8c) 


The additional diagrams arising from this expanded set of superheavy fermions are: 

Dll: 16, -1 -T26- 1 - V — 45\ -16,. I 61 -45 x - 1 '- 1 - lM - 1 - I 61 

Dll: I 61 -1-45^-1 - 1'-1-126-1 - I 61 , 

16i -1 -126 - 1 - V - 1 - 45^ -1 - 16! 

Dll: I 61 -1 - 45| -1 - 1' - 1 - 10' - 1 - 1' - 45 x -16 b 

16! —45 Y -V - 1 - 10' - 1 - V - 1 - 45^ -1 - I 61 

D 12 : I 61 -1 -T26 - 1 - V - 45 Y -16 2 

D 21 : 16 2 -45\- -V - 1 - 126 - 1 - I 61 

D 12 : I 61 -1-TM-45^-45 x -1 - 16 2 

D21: 16 2 -1 - 45 x -45^-126 - 1 - I 61 

D12: I 61 -1 - 45| -126 - 45 x -1 - 16 2 (3.7d) 
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D21: 16 2 —1 — 45x —126 — 45^ —1 — 16i 

D13: 16i -1 - 126 - 45* z -T* - 1 - 16 3 

D31: 16 3 -1 - 1'* - 45^-126 - 1 - 16i 

D13: 16i —1 — 45^ —126 — 1'* — 1 — 16 3 

D31: 16 3 -1 - 1'* -126 - 45^-1 - 16i 

D13: 16i -1 - 45| -1 - 10' - 45 x -16 3 

D31: 16 3 —45 y -10' - 1 - 45^ -1 - 16! 

Mil: 16! -45 x -T - 45^ -45^ -45 x -T - 126'- 

T - 45 y —45^ -45* z -T - 45 x -16, 

We thus have obtained the complete set of diagrams we shall consider for the evaluation 
of the mass matrices. Any additional diagrams for a given Mo irac or M R matrix element 
allowed by the U(1)f family symmetry are of higher-order and will lead to noticeably smaller 
contributions to that element than those arising from all the diagrams listed above. 

IV. EVALUATION OF CONTRIBUTIONS to the MASS MATRICES 

We now turn to the evaluation of the contributions to each matrix element at the SUSY 
GUT scale. The renormalizable 3-point couplings times VEVs for the 10(5), 10(5), 10 , (5 / ), 
126(5), 126(45) and 126 / (1) vertices contributing to M u and M N , M D and M E , M° and 
M e , M u and M N , M D and M E , and M R , respectively, are labeled 

9io v u, 9io v d, 9w v 'di 9i2eU! u , 9i26 w d, 9vm'W (4.1a) 

We shall assume the superheavy fermions all get massive at the same mass scale, so each 
1, 1', 45 Y or 45z vertex factor can be rescaled by the same propagator mass M according 
to 

X = 945 X U45 X /M, z = g 45z u 45z /M, s = giUi/M, s' = g v u v /M (4.1 b) 
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where we have introduced a convenient short-hand notation. In order to accommodate CP 
violation, as noted earlier after (3.5b) we introduce the four phases 

0 a: j 0z: 01: 01' (4.1c) 

As a result we are led to introduce 14 independent parameters in order to explain the 15 
quark and lepton masses and 8 quark and lepton mixing parameters. 

The contributions for each diagram then follow by moving along each fermion line and 
appending the above parameters together with the coupling coefhcents spelled out in Table I. 
Alternatively, one can use the detailed computational procedure of Mohapatra and Sakita 
[23] which makes explicit use of the SU( 5) decompositions of the S'O(IO) matrices and 
fields. We have used both procedures for a check in our calculations and both agree. In 
the expressions presented below, we have evaluated the Dirac (0 l) t C _1 ('0 c )l and Majorana 
('0 c )l T C'~ 1 ( , 0 c )l matrix elements . 

Leading-Order Dirac Matrix Diagrams of (3.7a): 

D33: 16 3 -10- 16 3 

Mg = Mg = 9w v ui Mg = Mg = 9io v d 

D23: 16 2 -10'- 16 3 

D32: 16 3 -10'- 16 2 

Mg = Mg = Mg = Mg = g' 10 v' d 

D22: 16 2 -126- 16 2 (4.2a) 

(-^22) -^ 22 ) = (1) ~ ^)9l26 w u: (Mg t Mg) = (1) ~‘i)9l2& w d 

D13: I 61 — 45 y -10' - 45 x -16 3 

D31: 16 3 — 45 y -10' - 45.y -I 61 

Mg = Mg = Mg = Mg = -3 g w v' d x 2 e 2i * x 
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D12: 16i -45 z -10'* - 45 z -16 2 
D21: 16 2 -45 z -10'* - 45 z -16i 

Mg = Mg = -4 g w v’ d z 2 e 2i *% Mg = Mg = 0 

Leading-Order Majorana Matrix Diagrams of (3.7b): 

M33: 16 3 -1 - TM' - 1 - 16 3 

Mg = gi26'w's 2 e 2l(t>1 

M23: 16 2 -1 - 45 A - -1' - TM' - 1 - 16 3 

M32: 16 3 -1 - TM' - V - 45 A - -1 - 16 2 

Mg = Mg = 5gi2&w'xs 2 s'e i ^ x+ ^ 1+ ^ 1 ^ 

M22: 16 2 -1 - 45 A - -1' - TM' - 1' - 45 x -1 - 16 2 

Mg = 25gi26'w' (xss') 2 e 2 ‘ l ^ x+ ^ 1+ ^ >1 '' > 

M13: 16i -45 a -1'-1-45 z -l'-lM' - 1 - 16 3 (4.2b) 

M31: 16 3 -1 - TM' - V - 45 z -1 - 1' - 45 x -16! 

Mg = Mg = 30 g u& w , xzs 2 s' 2 e i ^ + ^ +2 * 1+2 ^ 

M12: 16i —45 y -1' - 1 - 45 z -1' - TM' - 1' - 45 x -1 - 16 2 

M21: 16 2 -1 - 45 a - -1' - TM' - V - 45 z -1 - 1' - 45 x -16i 

Mg = Mg l = 150 gi26'w'x 2 zs 2 s' 3 e l ^ x+ ^ >z+2< ^ >1+3 ^ 

Mil: 16i -45 x -1' - 1 - 45 z -1' - TM' - 1' - 45 z -1 - V - 45 a -16i 

M/J = 9005'i 2 6'W / ( xzss 12 ) 2 e 2 *(<A+<A+</>i+2</>i/) 

Higher-Order Diagrams listed in (3.7c) from Minimal Set: 

16i —45 A -1' - 1 - 10' - 1 - 1' - 45 A -16i 

Mg = Mg = -Zgwv'gxss'Ye 2 ^^^^ 


Dll: 
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D22: 16 2 -45 z -10'* - V* - 16 2 , 16 2 — 1'* - 10'* - 45 z -16 2 

(4-2c) 

M23: 16 2 -45^ -1' - 1 - 45 z -1' - 126' - 1 - 16 3 

M32: 16 3 -1 - TM' - V - 45 z -1-1'- 45* -16 2 

= 30gi26>w , xzs 2 s' 2 e l( '-‘t >x+ '1 >z+2< t >1+2< t >1 ^ 

Higher-Order Diagrams of (3.7d) from the Expanded Set: 

Dll: 16i -1-126-1 - E-45 x -16i, 16i -45 x -1'- 1 - lM - 1 - 16i 

{M[[, Mg) = (2, -6)g 1 26W u xs 2 s'e i ^ +2 ^ +<t> ^ 

{Mg_, Mg) = (-2, 6)g 126 w d xs 2 s'e i ^ +2 ^ + ^ 

Dll: 16i -1 - 45^-1 - 1'- 1-126 - 1 - 16i, 

16i -1 - 126 - 1 - V - 1 - 45* -1 - 16! 

{Mg, Mg) = (3, —9)gi26W u zs 4 s , e i ^Mz+^4>i+4>i') 

{Mg, Mg) = (-3, 9 )g 126 w d zs 4 s'e i ^+Hi+M 

Dll: 16i -1 — 45^ —1 — 1' — 1 — 10' — 1 — 1' — 45* -16 u 

16! -45 x -1' - 1 - 10' - 1 - 1 ' - 1 - 45* z -1 - 16i 

{Mg, Mg) = (-7, -3)g w v' d xzs 4 s ,2 e i ^- <l ’z+Wi+ 2 M 

D12: 16i -1 -TM - 1 - V - 45* -16 2 

{Mg, Mg) = (1, —15) gi 26 w u xs 2 s'e^Mx+^4>i+4>i') 

Mg = Mg = ^3gi 2e w d xs 2 s'e^Mx+Wi+ti') 

D21: 16 2 —45* —V — 1 — 126 — 1 — 16i 

{Mg, Mg) = (1, 9)g 126 w u xs 2 s , e i ^ x+2 ^ + ^ 

{Mg, Mg) = (1, 9)g 12 6W d xs 2 s'e i (M*+Wi+<h’) 

D12: 16i —1 — TM — 45^ —45* —1 — 16 2 

{Mg, Mg) = (2, —90)gi 2e w u xzs 2 e i ( < t >x M*+ 2, t>i) 

Mg 2 = 12gi 26 w d xzs 2 e l ^ x ~ (t>z+2(t>1 \ Mg 2 = 0 
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D21: 16 2 -1 - 45 x -45^-126 - 1 - 16i 

(Mg, Mg) = (1, -27 )g 126 w u xzs 2 e i ^-^ +2 ^ 

{Mg, Mg) = (1, ~27)g V26 w d xzs 2 e i ^-^ +2 ^ ) 

D12: 16i -1 - 45^ -126 - 45 x -1 - 16 2 (4.2d) 

(Mg, Mg) = (1, 45 )g 12 6'w u xzs 2 e i ^- (t ’z+Wi) 

(Mg, Mg) = (-3, tyguewdxzs 2 ^^-^ 2 ^ 

D21: 16 2 -1 - 45 x -TM- 45^ -1 - 16i 

(Mg, Mg) = (2, 5A)g V26 w u xzs 2 e i ^- ,t, ,+^i) 

Mg = —4gi2 e w d xzs 2 e l< ' <t>x ~ <t>z+2<l>1 \ Mg = 0 

D13: 16i —1 — 126 — 45^ —l'* — 1 — 16 3 

(Mg, Mg) = (2, -18)g 126 wuzs 2 s' 

Mg = -Ag^w^s'e^Mz+WiMv), Mg = 0 

D31: 16 3 -1 - 1'* - 45^ -126 - 1 - 16i 

(Mg, Mg) = (1, 9)g 126 w u zs 2 s , e i ^+24>i-M 
(Mg, Mg) = (1, 9 )g 12 eW d zs 2 s'e ii -^+Wi-<fii') 

D13: 16i -1 - 45^ -126 - V* - 1 - 16 3 

(Mg, Mg) = ( 1 , 9 )g 1 2 Q WuZS 2 s'e i ^+ 2 ^-M 
(Mg, Mg) = (1, 9)gi26W d zs 2 s'e i (M*+ 2 4>i-<l>i>) 

D31: 16 3 -1 - 1'*-126-45^ -1 - 16i 

(Mg-y, Mg) = (2, -18)g 126 w u zs 2 s’ 

Mg = -4 g i 26 ^s 2 s / e i (-^ +2 ^-^'), Mg = 0 

16i -1 - 45^ -1 - 10' - 45 x —16 3 

Mg = Mg = —3gi 0 >v' d xzs 2 e i ^ x Mz+ 2 <t>i) 


D13: 
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D31: 


Mil: 


16:3 -45x -10' - 1 - 45* z -1 - 16! 

Mg = —4g 10 >v f d xzs 2 e i ^ x ~^ z+2 ^ >1 \ Mg = 0 

16! -45x -1' - 45^ -45* z -45 x -1' - 126'- 

1' - 45y -45| -45| -V - 45y -16i 

M/J = (900) 2 ^i26'W / (^s / ) 4 e 4i ^-^ + ^') 


An interesting observation which can be drawn from the Majorana contributions in (4.2b) 
is that the matrix in leading order has a geometrical texture as given in (2.4c) with 

Ml — 5 v /a!*M 3 r 3 (4.3) 

provided x ~ z. In fact, this observation served as an important guide in our construction of 
the Majorana neutrino matrix and suggested the relative roles played by the 45x and 45^ 
Higgs fields. 

V. QUANTITATIVE RESULTS for the 50(10) x U(1) F MODEL 

Finally we attempt to select a set of values for the 14 input parameters of (4.1a,b,c) 
which will accurately reproduce the input data in (2.1a,b) and (2.3a,b) used for our bottom- 
up approach. As noted earlier, the minimal set of superheavy fermions and their mirror 
partners found in (3.8b) yield unsatisfactory results: m u = rn Ve = 0, m e = 0.006 MeV 
and m„ = m„ T = 0.089 eV. The problem can be traced to the zero or tiny values of Dll. 
By expanding the set of superheavy fermions to include those in (3.8c), on the other hand, 
excellent results can be found as shown below. 

One particularly good numerical choice for the parameters at the SUSY GUT scale is 


given by 
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9w v u 

= 120.3, 

gwVd 

= 2.46, 

gwv’d 

= 0.078 GeV 

9l2§W u 

= 0.314, 

fh'XiWd 

= -0.037, 

912&w' 

= 0.8 x 10 16 i 

UAb x /M 

= 0.130, 

945z u 45z/M 

= 0.165, 

9iUi/M 

= 0.56, g\'U‘ 

(fix 

= 35°, 

4>z 

= 4>i = <t>v 

= -5° 



which reduces the number of independent parameters from 14 to 12. In fact, the only 
large phase angle is that for (j) x . As seen from (4.2a), this is in agreement with our earlier 
conclusion from the bottom-up phenomenological results [8] that essentially only the Dirac 
D13 and D31 matrix elements are complex. The mass matrices at the SUSY GUT scale are 
then numerically equal to 

/ -0.0010 -0.0001* 0.0053 + 0.0034i -0.0013 \ 


M u = | 0.0053 + 0.0034i 0.314 0 

-0.0013 0 120.3 

-0.0001 —0.0104 + 0.0004* -0.0029 - 0.0045*' 

-0.0077 + 0.0018* -0.036 0.078 

-0.0033 - 0.0048* 0.078 2.460 

’ 0.0030 + 0.0003* -0.079 - 0.051* 0.0038 \ 


V -0.0013 

/ - 0.0001 

M D = -0.0077 + 0.0018* 


M n = 0.048 + 0.031* 

V 0.0038 
/ 0.0004 

M e = 0.0060 + 0.0031* 


-0.942 0 

0 120.3 

-0.0020 - 0.0010* -0.0023 - 0.0045*' 
0.112 0.078 


(5.2a) 


(5.25) 


(5.2c) 


(5.2d) 


V-0.0009 - 0.0037i 0.078 2.460 ) 

( (-.069 + .640*) x 10 9 (-.141 - ,119*) x 10 11 (.108 + .019*) x 10 13 \ 

M R = (—.141 — .1197) x 10 11 (.461 + ,549i) x 10 12 (-.393 - .155*) x 10 14 (5.2e) 

V (.108 + .019*) x 10 13 (-.393 - ,155*) x 10 14 (.247 - .044*) x 10 16 


in units of GeV. By using the seesaw formula [16], we hnd for the light neutrino matrix at 
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the SUSY GUT scale 
M N *ff ~ -M n (M r )- 1 M nT 

(.027-,238i) x 10~ 3 (-.109 - .199*) x 10~ 2 (-.037 + .512*) x 10" 2 

(-.109 - .199*) x 10" 2 (-.232 - .088*) x 10” 1 (.258 + .435*) x 10” 1 

(—.037 + .512*) x 10“ 2 (.258 + .435*) x 10” 1 —.001 — .1127 

(5-2/) 

in electron-Volts. Again we emphasize the Dirac mass matrix elements appear in the form 
while the Majorana matrix elements refer to ('ip c )f L C~ 1 M(ip c )j] J with 
and (' 4> c )jL each a member of one of the three families of 16 s. Identical contributions also 
arise from the transposed Dirac matrices and the right-handed Majorana matrix. As such, 
the true Yukawa couplings Gy are just half the values of the gy’s appearing in (4.1a,b). 

The masses at the GUT scale can then be found by calculating the eigenvalues of the 
Hermitian product MAP in each case, while the mixing matrices Vckm and Vi ep t can be 
calculated with the projection operator technique of Jarlskog [24], After evolving these 
quantities to the low scale, we find in the quark sector 

m n (lGeV) = 5.0 (5.1) MeV, m d (lGeV) = 7.9 (8.9) MeV 

m c {m c ) = 1.27 (1.27) GeV, m s (lGeV) = 169 (175) MeV (5.3a) 

m t (m t ) = 150 (165) GeV, m b (m b ) = 4.09 (4.25) GeV 

where we have indicated the preferred values in parentheses. The mixing matrix is given by 

0.972 0.235 0.0037e- ,:i24 ° 

-0.235 0.971 0.041 

0.012 -0.039 0.999 

—0.003i -0.001* 

Note that V cb = 0.041 and \V ub /V cb \ = 0.090 with the CP-violating phase 5 = 124°, while 
m d/ m u — 1-59 and m s /rrid = 21.3, cf. [12, 13]. These results should be compared with our 
central starting input values given in (2.1a,b). 
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In the lepton sector we obtain 

m Ve = 0.10 (?) x 10 ~ 4 eV, m e = 0.43 (0.511) MeV 

m ViJ = 0.29 (0.25) x 10 " 2 eV, rri,, = 103 (105.5) MeV (5.4a) 

m Vr = 0.12 (0.10) eV, m T = 1.777 (1.777) GeV 

and 

/ 0.998 0.049 0.039e - * 121 

T r _ -0.036 0.875 0.483 

vlept — 

0.042 -0.482 0.875 

V —0.037i — 0 . 002 i 

which should be compared with the input values in (2.3a,b). The heavy Majorana neutrino 
masses are 



Mf = 0.63 x 10 9 GeV, Mf = 0.37 x 10 11 GeV, AFf = 0.25 x 10 16 GeV (5.4c) 


The neutrino masses and mixings are in the correct ranges to explain the nonadiabatic solar 
neutrino depletion with small mixing [5] and the atmospheric neutrino depletion with large 


mixing [ 6 ]: 


Sm\ 2 

= 8.5 x 10 " 6 eV 2 , 

sin 2 26*12 

= 0.0062 


= 1.4 x 10 ~ 2 eV 2 , 

sin 2 26*23 

= 0.71 


(5.5) 


For our analysis, the SUSY GUT scale at which the gauge and Yukawa couplings unify 
was chosen to be A = 1.2 x 10 16 GeV. From (3.4a) and (5.2a,b,c,d) we find that gw = 0.69. 
It is interesting to note that if we equate the S'O(10)-breaking and lepton number-breaking 
VEV, w', with A, we find gi 26 ' — 0.67 ~ gw Taking into account the remark following 
(5.2e), we note the true Yukawa couplings are Gho — — 0.33. If we further equate 

g L = g w — gi 2 &, and Ui = A for the [7(1 ^-breaking VEV, we End M = 1.5 x 10 16 GeV for 
the masses of the superheavy fermions which condense with their mirrors. These values are 
all very reasonable. 
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The 45\- and 45^ VEVs appear at nearly the same scale, 2.8 x 10 15 and 3.5 x 10 15 GeV 
respectively, if one assumes the same Yukawa coupling as above. On the other hand, if these 
VEVs appear at the unification scale A the corresponding Yukawa couplings are smaller than 
those found above. In either case, a consequence of their non-orthogonal breakings is that 
SU( 5) is broken down to SU( 3) c x SU(2)l x U(1)y at the scale in question. No further 
breaking is required until the electroweak scale and the SUSY-breaking scale are reached. 

VI. SUMMARY 

Our starting point for this research has been based on the results obtained from a bottom- 
up approach proposed previously by us to obtain mass matrices at the SUSY GUT scale 
based on a complete set of data inputted at the low scales. In particular we have used 
the known quark and charged lepton masses and CKM mixing matrix together with the 
neutrino masses and mixings based on particular neutrino scenarios. The masses and mixing 
matrices were evolved to the SUSY GUT scale where the mass matrices can be constructed 
by use of Sylvester’s theorem. By varying the bases and the signs of the mass eigenvalues, 
we looked for simple textures for the mass matrices such that each matrix element involved 
as few SO(IO) Higgs representations as possible. The neutrino scenario examined which 
appeared to yield the simplest structure involved the MSW nonadiabatic depletion of the 
solar electron-neutrinos together with the observed depletion of atmospheric muon-neutrinos 
by oscillations into tau-neutrinos. 

In this paper we have constructed an SO(10) x U(1)f model of the Yukawa interactions 
which neatly reproduces the desired SO(10) textures for the quark and lepton mass matrices 
for this preferred neutrino scenario. The observed features include the following: 

(i) The Abelian U(1)f family symmetry group singles out a rather simple set of tree 
diagrams which determines the texture of the generic Dirac and Majorana mass matrices, 
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while the 50(10) group relates corresponding matrix elements of the up, down, neutrino 
and charged lepton Dirac matrices to each other. 

(ii) The dominant second and third family Yukawa interactions are renormalizable and 
arise through couplings with Higgs in the 10, 10' and 126 representations of 50(10). The 
remaining Yukawa interactions are of higher order and require couplings of Higgs in the 
126 , 1, 1', 45 x and 45 z representations which acquire VEVs near the SUSY GUT scale. 

(iii) The Higgs which acquire high scale VEVs break the 50(10) xU(l)p symmetry down 
to the SU (3) c x SU (2 )l x U(1)y standard model symmetry in two stages through the SU (5) 
subgroup. 

(iv) Although this non-minimal supersymmetric model involves several Higgs represen¬ 
tations, the runnings of the Yukawa couplings from the GUT scale to the low-energy SUSY- 
breaking scale are controlled mainly by the contributions from the 10 , as in the minimal 
supersymmetric standard model. 

(v) The complete set of low scale VEVs which contribute to the fermion masses are 
10(5), 10(5), 10'(5'), 126(5) and 126(45) in the 5O(10)(5C/(5)) notation. These Higgs 
correspond to the minimum number required in 50(10) models which lead to the successful 
Georgi - Jarlskog relations [3]. Most of these models, however, do not include neutrino mass 
matrices. 

(vi) In terms of 12 input parameters, 15 masses (including the heavy Majorana masses) 
and 8 mixing parameters emerge. The Yukawa couplings and the Higgs VEVs are numeri¬ 
cally feasible and successfully correlate all the quark and lepton masses and mixings in the 
scenario which incorporates the nonadiabatic solar neutrino and atmospheric neutrino de¬ 
pletion effects. 

(vii) The right-handed Majorana neutrino matrix has a nearly geometrical texture lead¬ 
ing to heavy Majorana neutrino masses spread over seven orders of magnitude as given in 
(5.4c). In fact, it is the highly geometrical structure of the Majorana matrix which accounts 
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for the nearly maximal mixing of the ry and u T , rather than sizable mixing in the Dirac 
sector [25]. 

With the model as presented, the U(l)p current is anomalous, since the [f/(l),p] 3 triangle 
anomaly does not vanish, ft is possible to cancel this anomaly, however, by the addition 
of two 50(10) singlet neutral fermions, n i and ( u c )l , both with U(1)f charges of -12. By 
introducing another Higgs singlet representation which develops a GUT scale VEV, one can 
arrange that one of the new neutrinos remains massless while the other becomes superheavy. 
Alternatively, it is possible to cancel such an anomaly through the Green-Schwarz mechanism 
[22] provided the model can be derived from string theory. 

Studies are underway to examine what effects small mixings of such a light sterile neu¬ 
trino with the three families of light neutrinos will have on the neutrino spectrum and will 
be reported elsewhere. Work is also underway to construct a superpotential for the model 
presented here. 
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SU(5) 

Assignments 

VEV Directions 

45 x 45 y 45 z 

Flipped SU(5) 

Assignments 

u, d 

1 

1 

1 

d, u 

u c 

1 

- 4 

2 

d c 

d c 

- 3 

2 

- 4 

u c 

V, i 

- 3 

- 3 

- 3 

£, v 

v c 

5 

0 

6 

e c 

e c 

1 

6 

0 

v c 


Table I. Couplings of the 45 VEVs to states in the 16. 


Fig. 1. Tree-level diagrams for the (a) renormalizable and (b) leading-order nonrenormal- 
izable contributions to the generic Dirac mass matrix and for the (c) 33 element of the 
Majorana mass matrix. 
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